A simplified approach for Chaplygin-type cosmologies 
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A new class of accelerating cosmological models driven by a one-parameter version of the general 
Chaplygin-type equation of state is proposed. The simplified version is naturally obtained from 
causality considerations with basis on the adiabatic sound speed 115 plus the observed accelerating 
stage of the universe. We show that very stringent constraints on the unique free parameter a 
describing the simplified Chaplygin model can be obtained from a joint analysis involving the latest 
SNe type la data and the recent Sloan Digital Sky Survey measurement of baryon acoustic oscilla- 
tions (BAO). In our analysis we have considered separately the SNe type la gold sample measured 
by Riess et al. (2004) and the Supernova Legacy Survey (SNLS) from Astier et al. (2006). At 95.4% 



(c.L), we find for BAO + gold sample, 0.91 < a < 1.0 and Am 
analysis provides 0.94 < a < 1.0 and Qm 



0.28 



+0.043 



while BAO + SNLS 



n 97 +0.048 
u - z ' -0.045- 



PACS numbers: 98.80.Es; 95.35.+d; 98.62.Sb 



I. INTRODUCTION 

The impressive convergence of recent observational 
facts along with some apparently successful theoretical 
predictions seem to indicate that the simple approach 
provided by the standard cold dark matter (CDM) model 
is insufficient to describe the present stage of our uni- 
verse. From these results, the most plausible picture for 
our world seems to be a nearly flat scenario dominated 
basically by CDM and an exotic component endowed 
with large negative pressure, usually named dark energy. 
Despite the good observational indications for the exis- 
tence of these two components, their physical properties 
constitute a completely open question at present, which 
gives rise to the so-called dark matter and dark energy 
problems (see pj for a recent review on this topic). 

Among the many candidates for the dark energy com- 
ponent, a very interesting one was suggested by Kamen- 
shchik et al. 0| and developed by Bilic et al. Q and 
Bento et al. 4,] . Such an exotic fluid, named generalized 
Chaplygin gas (C-gas), can be macroscopically charac- 
terized by the equation of state (EoS) 



PC 



-A/p a c , 



(1) 



where a — 1 and A is a positive constant related to 
the present-day Chaplygin adiabatic sound speed, v 2 s = 
aA/p^ a (pc stands for the current C-gas density). In 
actual fact, the above equation for a =/= 1 constitutes 
a generalization of the original C-gas EoS proposed by 
Bento et al. in Ref. 0. 
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In the last few years, the possibility of describing the 
unknown dark energy component using the C-gas-type 
EoS above has provoked a considerable debate in the lit- 
erature. Theoretical connections between the C-gas and 
string theory, supersymmetric generalizations |3, |6| , self- 
interacting [3 , and even a tachyonic fluid representation 
Q has also been investigated. Another interesting fea- 
ture of the above EoS comes from the fact that the C-gas 
becomes pressureless at high redshifts, which suggests a 
possible unification scheme for the cosmological "dark 
sector" , an interesting idea which has been considered in 
different contexts [jj. 

Observational aspects of the above C-gas scenarios 
have also been largely investigated in the literature. Cos- 
mological tests involving type la supernovae (SNe la) 
data |lOl ITU , the shape of the matter power s pec trum 
[l2| , statistical properties of gravitational lenses |13| , the 
age of the Universe [DJ, cosmic microwave background 
(CMB) measurements |l5l llH Il7t EH , galaxy clusters X- 
ray [l!j , and gamma-ray bursts data [2(| have been dis- 
cussed. In general, to perform such analyses, besides the 
present value of the C-gas density parameter (^c)i the 
above barotropic EoS implies that one needs to constrain 
two additional free parameters, namely, A and a. There- 
fore, in the context of the Friedman-Robertson- Walker 
(FRW) cosmologies with CDM plus a C-gas, there are at 
least 4 parameters to be constrained by the data. Actu- 
ally, this number can be reduced to 3 if one assumes a 
flat geometry, i.e., £Im = 1 — or if a unified dark mat- 
ter/energy picture involving only the C-gas and baryons 
is assumed from the very beginning (in this case, the 
baryonic density (f2f,) may be fixed a priori by using, 
for instance, nucleosynthesis or the recent Cosmic 
Microwave Background (CMB) observations [2j|). How- 
ever, even in this latter cases, there are so many parame- 
ters to be constrained by the data, that a high degree of 
degeneracy on the parametric space becomes inevitable. 



2 



Many generalizations of the original C-gas [24L 125. 26, 
|23,|1^|, or even of its extended version |2jJ have appeared 
in literature. In these cases, the number of free parame- 
ters is usually increased, and, as consequence, the models 
become mathematically richer although much less predic- 
tive from a physical viewpoint. In this work by following 
the opposite direction, we propose a simplified version 
for the generalized C-gas-type EoS which diminishes one 
of its free parameters. By an additional physical con- 
dition, the allowed range of the remaining parameter is 
also restricted a priori, which makes not only the rele- 
vant parametric space bi-dimensional but also (and more 
important) the model more easily discarded or confirmed 
by the present set of observations since the range of its 
free parameter is physically limited from causality con- 
siderations. We test the viability of this simplified C-gas 
approach by discussing the constraints imposed from cur- 
rent SNe la observations and Large Scale Structure (LSS) 
data. 



II. A SIMPLIFIED C-GAS SCENARIO 

Let us consider a homogeneous and isotropic Universe 
whose energy components are cold dark matter plus the 
generalized C-gas fluid. Since both components are sep- 
arately conserved, by inserting Eq. (1) into the energy 
conservation law pc = —3H(pc +Pc), one obtains the 
following expression for the density of the C-gas 0, 0, 0] 



PC = PC 



A s + {1-A s )a 



3(l+a) 



(2) 



where a(t) is the cosmological scale factor and A s = 
Aj ' p X Q~ a is a convenient dimensionless constant (as usual, 
the subscript "0" denotes present-day quantities) . As one 
may check, the above C-gas evolving in the FRW met- 
ric can be modeled as a quintessence, that is, a scalar 
field model described by an ordinary Lagrangian density, 



V ((/>), with the following potential 



V{4>) = ^Pc As +1 {[coshV6nm pl (a + l)(j))] — 

2a 



[cosh v Q-Km p i(a + 1) 



(3) 



where m p i is the Planck mass. 

In a flat geometry, the Friedmann equation for a con- 
served C-gas plus cold dark matter is given by 



n 



n M a~ 3 + fl c [A s + (1 - A s )a 



-3(a+l) 



(4) 



where U = H(a) 2 /H$. Note that besides the Hubble 
parameter Hq we still have 3 additional parameters in 
this case (a, A s , Qm), even using the flat condition Qc — 
1 — Qm- Therefore, an interesting question to be answered 
at this point is how to reduce the C-gas parameters based 
on reasonable physical constraints? 

In order to answer the above question, we first no- 
tice that the dimensionless constant A s appearing in the 



above expressions encodes the basic information coming 
from the original parameter A [see Eq. (1)] . On the other 
hand, the Chaplygin adiabatic sound speed reads 



dp 
dp 



v 2 = £ = aA/p l c+ « 



(5) 



which must positive definite for a well-behaved gas (zero 
in the limit case of dust). Note also that the present day 
Chaplygin adiabatic sound speed is v^ — aA/ p c i+ a or, 
equivalently, 



v 2 so = aA/p 1 ^ = aA s . 



(6) 



Therefore, if the A s parameter is a function of a, the 
number of free parameters is naturally reduced, and, as 
an extra bonus, the positiviness of v 2 at any time, as well 
as its thermodynamic stability, is naturally guaranteed. 
Among many possible relations (e.g., A s — a n ), clearly 
the simplest choice is A s = a (n = 1). In this case, 
v 1o = a2 > or more generally, v 2 = a 2 (pco/ p) a ■ Note 
also that, since the light speed is a natural cutoff for 
the sound propagation, it follows that v so — \a\ < 1, 
thereby restricting a to the interval [-1,1]. An additional 
constraint can still be imposed to this parameter. In fact, 
with A s — a, the simplified C-gas EoS (1) becomes 



Pc 



-apco 



PCo 

Pc 



(7) 



so that a negative pressure is obtained only for positive 
values of a. In other words, this accounts to saying that 
the combined requirements from causality along with the 
observed accelerating stage of the Universe limit natu- 
rally the parameter a to the interval < a < 1. 

Note that the the simplified Chaplygin gas above (from 
now on SC-gas) preserves the unifying character of the 
original C-gas, i.e., it behaves as a pressureless fluid (non- 
relativistic matter) at high-z while, at late times, it ap- 
proaches the quintessence behavior, which now is fully 
characterized by the a parameter (for a unified dark 
matter / dark energy description of the above scenario, 
see |21|). However, note also that, even in this limit- 
ing case, the sound speed is positive. In other words, 
the Universe evolution resembles the one driven by a 
quintessence component but the thermodynamic behav- 
ior does not present the pathologies of such scenarios. 

In this simplified approach, Eq.QJ is rewritten as 



n 



n M a- 3 + n c [a + (1 - a)a- 3(a+1) ] 



(8) 



so that the parameter a is actually the unique unknown 
constant related to this SC-gas model. In what follows, 
we confront this simplified approach with the most recent 
SNe la and Large Scale Structure (LSS) data. 
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FIG. 1: Supernova results. Panel (a) displays the residual magnitude with respect to an empty (Qt = 0) universe for the 
HZST (filled circles) and SNLS (open circles) samples, respectively. The solid curves are the predictions for the SG-gas-type 
models characterized by (0m, a) • F° r comparison, we also display the predictions of the open cold dark matter (OCDM) and 
Einstein-de Sitter Universes. Panels (b) and (c) show 68.3%, 95.5% and 99.73% confidence contours on the space parameter 
(f2M,a) from the HZST and SNLS supernova data, respectively. The horizontal dotted and dashed lines correspond to the 
constraints arising from the SDSS baryon acoustic oscillations detection. For SNe la data alone, if a is greater than 0.91, values 
of ilM smaller than 0.15 are ruled out by the two samples. Constraints from BAO contribute to increase the allowed values of 

Om- 



ul OBSERVATIONAL CONSTRAINTS 
A. SNe la 

Let us first investigate the bounds arising from SNe 
la observations on the SC-gas scnario described above. 
To this end we use the most recent SNe la observa- 
tions, namely, the High-Z SN Search (HZS) Team [U 
and the Supernova Legacy Survey (SNLS) Collaboration 
data H3. 

The so-called gold sample from the HZS team is a se- 
lection of 157 SNe la events distributed over the red- 
shift interval 0.01 < z < 1.7, and constitutes the com- 
pilation of the best observations made so far by them 
and by the Supernova Cosmology Project plus 16 new 
events observed by Hubble Space Telescope (HST). The 
current data from SNLS collaboration correspond to the 
first year results of its planned five year survey. The to- 
tal sample includes 71 high-z SNe la in the redshift range 
0-2 < z < 1 plus 44 low-z SNe la. This data set is ar- 
guably (due to multi-band, rolling search technique and 
careful calibration) the best high-z SNe la compilation 
to date, as indicated by the very tight scatter around the 
best fit in the Hubble diagram and a careful estimate of 
systematic uncertainties. Another important aspect to 
be emphasized on the SNLS data is that they seem to be 
in a better agreement with WMAP results than the gold 
sample (see, e.g., for a discussion). The two SNe 



la samples are illustrated on a residual Hubble Diagram 
with respect to the empty universe model (f2x = 0) in 
Fig. la. 

The predicted distance modulus for a supernova at red- 
shift z, given a set of parameters p, is 

MpOIp) = m - M = 51ogd L + 25, (9) 

where m and M are, respectively, the apparent and ab- 
solute magnitudes, the complete set of parameters is 
p = (H a , Om, ol) and g?l stands for the luminosity dis- 
tance (in units of megaparsecs), 

with x' = (1 + z) -1 being a convenient integration vari- 
able and Tt(x; p) the expression given by Eq. ©. 

We estimated the best fit to the set of parameters p 
by using a \ 2 statistics 

2 _ ^ [/4(^|p) -MoOIp)] 2 , s 
x — 2 ' \ LL ) 

with the parameters f^M and a spanning the interval [0,1] 
in steps of 0.01. In the above expression, N = 157 and 
115 for gold and SNLS samples, respectively, ii l p (z\p) is 
given by Eq. 10, ^(z|p) is the extinction corrected dis- 
tance modulus for a given SNe la at Zi, and o~i is the 
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FIG. 2: Confidence contours on the space parameter (Qm, a) from a joint analysis involving SNe la and the SDSS baryon acoustic 
oscillations. As in Figure 2, panel (a) display the results for the gold sample while Panel (b) for the SNLS collaboration. The 
important aspect here is that the parameter a of the simplified model is extremely restricted by both set of data. 



uncertainty in the individual distance moduli. In our 
analysis, Hq is considered a nuisance parameter so that 
we marginalize over it. 

In Figures (lb) and (lc) we show the results of our 
statistical analysis. Contours of constant likelihood 
(99.73%, 95.4% and 68.3%) are shown in the paramet- 
ric space a — Qm- Panel (lb) displays the results for the 
HZS gold sample. Compared to Fig. 4 of Ref. ^lj, the 
parameter a is now considerably more restricted than in 
the standard C-gas approach. In particular, note that 
for any value of the matter density parameter, mod- 
els with a < 0.63 are ruled out at 99.73% level. The 
best-fit model for this analysis occurs for Q m = 0.0 and 
a = 0.79 with Xmin/^ = 1-13 (y = degrees of freedom). 
At 95.4% c.l. we also find Q M < 0.36 and 0.71 < a < 1.0. 
Panel (lc) shows a similar analysis for the SNLS data. 
The best-fit parameters in this case are Qm — 0.2 and 
q = 0.96 with xLin/ 1 ' — 1-0- Note that, when com- 
pared with recent dynamical estimates of Qm, this lat- 
ter value for the matter density parameter seems to be 
more realistic than the one provided by the gold sample 
analysis. The SNLS sample also imply Qm < 0.34 and 
0.75 < a < 1.0 at 95.4% (c.l.). 



B. SNe la + LSS analysis 

The recent detection of a peak in the lar ge s cale cor- 
relation function at lOO/i" 1 Mpc separation [33] provide 
not only a remarkable confirmation of the big bang cos- 
mology but also a kind of "ruler" with which cosmolog- 
ical scenarios can be tested. The peak detected (from 
a sample of 46748 luminous red galaxies selected from 
the SDSS Main Sample) is predicted to arise precisely at 
the measured scale and is basically due to baryon acous- 



tic oscillations (BAO) in the primordial baryon-photon 
plasma prior to recombination. Here, this measurement 
is characterized by 



A = 



Q 



1/2 
M 



1/3 



2/3 



0.469 ±0.017, (12) 



where z* = 0.35 is the redshift at which the acoustic 
scale has been measured, and r(z„) is the dimensionless 
comoving distance to z*. 

The dotted lines in Figs, (lb) and (lc) represent the 
constraints from SDSS BAO measurements on the pa- 
rameter space Qm — ct. Note that they are approxi- 
mately orthogonal to those arising from SNe la data, 
which indicates that possible degeneracies in the Qm — ct 
plane may be broken from a joint analysis involving 
these observational data sets. This is exactly what we 
show in Panels (2a) and (2b) for the BAO+gold and 
BAO+SNLS samples, respectively. Note that the avail- 
able parametric plane in both cases is considerably re- 
duced relative to the former analyses (Figs, lb and lc). 
Note also that, although compatible with the data, the 
region a > 1 (forbidden from thermodynamic stability 
and causality considerations) should be disregarded from 
the analysis since these arguments lead to the physi- 
cal bound < a < 1. For the BAO+gold sample we 
find Q M = 0.28±g g^| and a > 0.916 (with the best fit 
a = 0.98) at 95.4% (c.l.) while for the BAO+SNLS sam- 
ple the best-fit model happens at Qm = 0.27 and a = 1.0. 
This latter best-fit scenario corresponds to an accelerat- 
ing universe with go — —0.5, a total age of the Universe 
of t a ~ 10.2/i -1 Gyr, and a D/A redshift transition (from 
deceleration to acceleration) z d /a — 0.75. At 95.4% c.l., 
the BAO+SNLS analysis also provides 0.94 < a < 1.0 
and Q M = 0.27±g;gg. 



5 



IV. DISCUSSION AND CONCLUSIONS 

As widely known, there are many theoretical ap- 
proaches for describing the exotic dark energy component 
accelerating the Universe. However, until the present, 
the available battery of cosmological tests was not ca- 
pable to decide which is the best theoretical representa- 
tion. We have argued here that one of such candidates, 
the so-called Chaplygin type gas (whose equation of state 
depends on two parameters A s and a), may have a very 
simplified description. We postulate that A s is a function 
of a and for simplicity we have taken A s = a. Thus, sim- 
ilarly to the concordance model (ACDM), the resulting 
flat cosmology is completely described only by a pair of 
parameters (a, Qu)- This SC-gas cosmology mimics the 
dynamics of the X-matter models with an extra bonus, 
namely: the fluid stability and other thermodynamic fea- 
tures are guaranteed from the very beginning. 

By considering this particular parameterization we 
have investigated constraints on the a parameter from 
the most recent SNe la (gold and SNLS samples) and LSS 



data. We have found that the limits arising from this par- 
ticular combination of the data are much more restrictive 
on this simplified approach than on the generalized C-gas 
version. In particular, for the the BAO+SNLS combina- 
tion we found 0.94 < a < 1.0 and fl M = 0.27^;^ (at 
95% c.L), which is in agreement with recent estimates of 
the clustered matter. Naturally, it should be interesting 
to investigate whether current CMB data and other inde- 
pendent observations can or cannot discard the simplified 
scenario proposed here. 

Acknowledgments 

This work is partially supported by the Conselho 
Nacional de Desenvolvimcnto Cienti'fico e Tecnologico 
(CNPq - Brazil). JASL and JVC thanks FAPESP 
(Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
- Brazil). JSA is also supported by Fundagao de Amparo 
a Pesquisa do Estado do Rio de Janeiro (FAPERJ) No. 
E-26/171. 251/2004. 



[1] T. Padmanabhan, Phys. Rep. 380, 235 (2003); V. Sahni 
and A. A. Starobinsky, Int. J. Mod. Phys. D9, 373 (2000); 
P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 

(2003) ; J. A. S. L ima, Braz. J. Phy s. 34, 194 (2004); E. 
J. Copeland et al, hep-th/0603057 

[2] A. Kamenshchik, U. Moschella and V. Pasquier, Phys. 

Lett. B 511, 265 (2001). 
[3] N. Bilic, G. B. Tupper and R. D. Viollier, Phys. Lett. B 

535, 17 (2002). 
[4] M. C. Bento, O Bertolami and A. A. Sen, Phys. Rev. 

D66, 043507 (2002). 
[5] M. Bordemann and J. Hoppe, Phys. Lett. B317, 315 

(1993); See also, J. Hoppe, hep-th/9311059. 
[6] R. Jackiw, physics/0010042. 

[7] H. Zhang and Z.-H. Zhu, Phys. Rev. D73, 043518 (2006). 
[8] See for instance, H. B. Benaoum, hep-th/0205140, and 

articles quoted there. 
[9] T. Matos and L. A. Urena-Lopez, Class. Quantum Grav. 

17, L75 (2000); Phys. Rev. D 63, 063506 (2001); A. 

Davidson, D. Karasik and Y. Lederer, gr-qc/0111107 C. 

Watterich, Phys. Rev. D 65, 123512 (2002); S. Kasuya, 

Phys. Lett. B 515, 121 (2001); T. Padmanabhan and T. 

R. Choudhury, Phys. Rev. D 66, 081301 (2002); M. Nov- 

ello, M. Makler , L.S. Werneck and C.A. Romero, Phys. 

Rev. D71, 043515 (2005). 
[10] M. Makler, S. Q. de Oliveira and I. Waga, Phys. Lett. B 

68, 123521 (2003); R. Colistete Jr., J. C. Fabris, S. V. B. 

Gongalves, P. E. de Souza, Int. J. Mod. Phys. D13, 669 

(2004) ; O. Bertolami, A. A. Sen, S. Sen and P. T. Silva, 
Mon. Not. Roy. Astro. Soc. 353, 329 (2004). 

[11] J. S. alcaniz an d J. A. S . Lima, Astrophys. J. 618, 16 

(2005) . astro- ph/0308465 

[12] P. P. Avelino, L. M. G. Bega, J. P. M. de Carvalho, C. 
J. A. P. Martins and P. Pinto, Phys. Rev. D 67, 023511 
(2003); L. M. G. Bega, P. P. Avelino, J. P. M. de Carvalho 
and C. J. A. P. Martins, Phys. Rev. D 67, 101301 (R) 



(2003) ; R. R. R. Reis, M. Makler and I. Waga, Phys. 
Rev. D69, 101301 (2004). 

[13] A. Dev, J. S. Alcani z and D. Ja in, Phys. Rev D 67, 
023515 (2003). astro-ph/0209379 P. T. Silva and O. 
Bertolami, Astrophys. J. 599, 829 (2003); A. Dev , D. 
Jain and J. S. Alcaniz, Astron. Astrophys. 417, 847 

(2004) . e-Print Archive: astro-ph/0311056 

[14] J. S. Alcaniz, D. Jain and A. D ev, Phys. Rev. D 67, 

043514 (2003). astro-ph/0210476 
[15] M. C. Bento, O. Bertolami and A. A. Sen, Phys. Rev D 

67, 063003 (2003) 
[16] M. C. Bento, O. Bertolami and A. A. Sen, Phys. Lett. B 

575, 172 (2003); Gen. Rel. Grav. 35, 2063 (2003). 
[17] D. Carturan and F. Finelli, Phys. Rev. D 68, 103501 

(2003). 

[18] L. Amendola, F. Finelli, C. Burigana and D. Carturan, 

JCAP 0307, 005 (2003). 
[19] J. V. Cunha, J. S. Alcaniz and J. A. S. L ima, Phys. 

Rev. D69, 083501 (2004). astro-ph/0306319 Z.-H. Zhu, 

Astron. Astrophys. 423, 421 (2004). 
[20] O. Bertolami and P.T. Silva, Mon. Not. Roy. Astron. Soc. 

365, 1149 (2006). 
[21] J. V. Cunha, J. S. Alcaniz and J. A. S. Lima, Submitted 

for Publication. (2006). 
[22] G. Steigman, Int. J. Mod. Phys. E 15, 1 (2006); J. 

P. Kneller and G. Steigman, Phys. Rev. D67, 063501 

(2003). 

[23] D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003); 

D. N. Spergel et al, 2006. astro-ph/0603449. 
[24] X. Zhang, F.-Q. Wu and J. Zhang, J. Cosmol. Astropart. 

Phys. 0601, 003 (2006). 
[25] L. P. Chimento and R. Lazkoz, Phys. Lett. B, 615, 146 

(2005) . 

[26] M. K. Mak and T. Harko, Phys. Rev. D71, 104022 
(2005). 

[27] Z.-K. Guo and Y.-Z. Zhang, astro-ph/0506091 



6 



[28] W. Wang et al., Mod. Phys. Lett. A20, 1443 (2005). 
[29] A. A. Sen and R. J. Scherrer, Phys. Rev. D72, 063511 
(2005). 

[30] A. G. Riess et al, Astrophys. J. 607, 665 (2004); 
[31] P. Astier et al, Astron. Astrophys. 447, 31 (2006) 



[32] H. K. Jassal, J . S. Bagla, and T. Padmanabhan, 

astro-ph/0601389 
[33] D. J. Eisenstein et al., Astrophys. J. 633, 560 (2005). 



